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Summary 



We present a general approach to the computational design of nanostructured 
chemical sensors. The scheme is based on identification and calculation of mi- 
croscopic descriptors (design parameters) which are used as input to a thermo- 
dynamic model to obtain the relevant macroscopic properties. In particular, we 
consider the functionalization of a (6,6) metallic armchair single-walled carbon 
nanotube (SWNT) by nine different 3d transition metal (TM) atoms occupying 
three types of vacancies. For six gas molecules (N2, O2, H2O, CO, NH3, H2S) we 
calculate the binding energy and change in conductance due to adsorption on each 
of the 27 TM sites. For a given type of TM functionalization, this allows us to ob- 
tain the equilibrium coverage and change in conductance as a function of the par- 
tial pressure of the "target" molecule in a background of atmospheric air. Specifi- 
cally, we show how Ni and Cu doped metallic (6,6) SWNTs may work as effective 
multifunctional sensors for both CO and NH3. In this way, the scheme presented 
allows one to obtain macroscopic device characteristics and performance data for 
nanoscale (in this case SWNT) based devices. 
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1 Introduction 



Detecting specific chemical species at small concentrations is of fundamental im- 
portance for many industrial and scientific processes, medical applications, and 
environmental monitoring. Nanostructured materials are ideally suited for sensor 
applications because of their large surface to volume ratio, making them sensitive 
to the adsorption of individual molecules. 

At a general level, any nanosensing system consists of the following four main 
components: 

1. a "target molecule" to be detected, 

2. an "active site" where the target molecule may adsorb on the sensor, 

3. a "sensing property" which changes upon adsorption of the target molecule, 

4. a "background" of adsorbing molecules which make up the background sig- 
nal. 

The active site must be designed so that adsorption of the target molecule in the 
presence of the background is sufficient to change the sensing property. These 
four main components of a chemical sensor are shown schematically in Figure [TJ 
for the case of a transition metal (TM) doped (6,6) metallic armchair single- waUed 
carbon nanotube (SWNT) measuring CO under atmospheric conditions. 

Performing a screening study of possible SWNT dopings experimentally would 
be quite demanding, especially considering the great variety of potential active 
sites and target molecules one may wish to detect. On the other hand, compu- 
tational screening studies based on density functional theory (DFT) calculations 
have recently been shown to be quite effective at shortening lists of potential can- 
didates, which may then be studied experimentally. 

The energetics obtained from DFT calculations using a generalized gradient 
approximation (GGA) for the exchange and correlation (xc)-functional are suffi- 
ciently accurate to provide a quantitative description of the stability of an active 
site, and the adsorption energies for both the target molecule and the background. 
From this data we may use kinetic modeling to estimate the equilibrium coverage 
of active sites by target molecules as a function of their pressure in the presence 
of the background. This allows us to quickly screen a wide variety of active sites, 
and rule out those which are completely oxidized in air, or too unstable at normal 
temperatures. 
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sensing property 



Figure 1: Schematic of a chemical sensor consisting of active sites (metal dopants 
in a (6,6) SWNT), a target molecule (CO), a background (atmospheric air), and a 
sensing property (resistance) (Mowbray et al. 2010). 

At the same time, changes in a nanosensor's resistance are often used to de- 
tect adsorption of a target molecule. Trends in the resistance may be reasonably 
described using the non-equilibrium Green's function (NEGF) methodology. In 
this way we may estimate the resistance per active site in the background, and 
hence determine the change in resistance of the active site as a function of target 
molecule concentration. 

Using resistance as a sensing property has two main advantages. First, it is a 
non-intrusive measure, which should not significantly influence the adsorption of 
the target molecules. Second, small changes in the coverage of target molecules 
(~ 0.1%) often yield large changes in resistance (~ IQ/site). 

Overall, this procedure fits together like four "puzzle pieces", DFT calcula- 
tions, kinetic modeling, NEGF results and circuit theory, to calculate the average 
change in resistance AR as a function of a target molecule's concentration C[X]. 
This methodology for modeling the chemical sensor is depicted schematically in 
Figure [2j 

From DFT we obtain the formation energies fiform of an active site (a 3d TM 
occupied vacancy (TM@VC)), the adsorption energies Eads for species X on the 
active site (X-l-TM@VC), and the Hamiltonian ^ and overlap matrices S for the 
scattering region with X adsorbed (X-l-TM@VC) and the leads (a pristine (6,6) 
SWNT (NT)). The energetics from DFT are then inputted into a kinetic model 
for the active site in thermodynamic equilibrium with concentrations Cq, and gas 
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DFT 



Kinetics 
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Figure 2: Schematic of the overall procedure for modeling a nanosensor, fitting 
together like four pieces of a puzzle, density functional theory (DFT) calcula- 
tions, kinetic modeling, non-equilibrium Green's function (NEGF) calculations, 
and circuit theory, as described in the text. 

phase entropies 5 gas, for each species taken from experiment (Lide 2006-2007). 
From this we obtain the fractional coverage per active site (TM@VC) as a 
function of the target molecule's concentration C[X]. 

At the same time, the Hamiltonian ^ and overlap 5 matrices obtained from 
DFT are used within the NEGF methodology to obtain the density of states (DOS) 
and transmission probability T for each species X on an active site (X-l-TM@ VC), 
as a function of energy s. Finally the coverage and transmission probabilities 
are used within simple circuit theory to estimate the scattering resistance for 
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species X on an active site (X+TM@VC), and hence derive the overall change 
in resistance AR for a particular active site (TM@VC) as a function of the target 
molecule's concentration C[X]. 

As a specific example of a chemical nanosensor we will consider the function- 
alization of a (6,6) SWNT by nine different 3d TM atoms occupying three types 
of vacancies. For six gas molecules (N2, O2, H2O, CO, NH3, H2S) we calculate 
the binding energy and change in conductance due to adsorption on each of the 
27 TM sites. For a given type of TM functionalization, this allows us to obtain 
the equilibrium coverage and change in conductance as a function of the partial 
pressure of the "target" molecule in a background of atmospheric air, as shown 
schematically in Figure [Tj and described in (Garcfa-Lastra et al. 2010, Mowbray 
et al. 2010). 

We will begin by providing a brief overview of the properties of our candi- 
date sensing materials, specifically 3d TM doped (6,6) SWNTs. We will then 
discuss in detail in separate sections the four components of our method for mod- 
eling nanosensors: DFT calculations, kinetic modeling, NEGF results, and circuit 
theory, along with the results obtained for the TM functionalized SWNTs. 



2 TM Doped SWNTs as Nanosensors 

To understand what a SWNT is (lijima 1991, Sanchez-Portal et al. 1999, Dres- 
selhaus et al. 2001, Harris 1999, Chalier et al. 2007), one may imagine a sheet of 
graphene being rolled into a diminutive cylinder, which is very long (^ ~ 1 cm) 
and narrow (J ~ 1 nm). The properties of the SWNT strongly depend on the 
way the graphene sheet is wrapped. To represent the way the graphene sheet is 
wrapped, two indices are used: (n, m), where n and m are integers. 

These two indices represent the number of unit cell vectors along the two 
principle directions in the crystal lattice of graphene, ai and 82. For a (n,m) 
SWNT, the circumference C is then given by 

C = llnai +ma2l|. (1) 

The nanotube is metallic if n- m = 3p, where p is an integer, and semiconducting 
otherwise. If n = m the tube is classified as an armchair SWNT, if m = the 
tube is called a zigzag SWNT, and otherwise the tube is known as chiral. Based 
on simple geometrical arguments, the diameter of a SWNT may be approximated 
using the relation 

d ~ — Vn^ -I- nm + nP-, (2) 
n 
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where a w 2.46 A. The strength, electronic structure, and high surface to volume 
ratio of SWNTs make them excellent materials for nanodevices. 

SWNTs work remarkably well as detectors for small gas molecules, as demon- 
strated for both individual SWNTs (Kong et al. 2000, Collins et al. 2000, Hierold 
2008, Villalpando-Paez et al. 2004, Rocha et al. 2008, Brahim et al. 2008) and 
SWNT networks (Morgan et al. 2008, Mowbray et al. 2009, Goldoni et al. 2004). 
However, such sensitivity is most likely attributable to structural defects, vacan- 
cies, and junctions in the SWNT networks, as pristine SWNTs are inherently inert 
(Mowbray et al. 2009). Thus, to control a SWNT's performance as a chemical 
sensor, one must control the location, type, and number of these structural de- 
fects. For this reason, the controlled doping and functionalization of SWNTs has 
become an area of increasing interest (Ayala et al. 2010), especially for substitu- 
tional doping with either N (Ayala et al. 2007) or B (Ayala et al. 2008). 

Previous studies have shown that SWNTs are highly sensitive to most molecules 
upon functionalization (Fagan et al. 2003, Yagi et al. 2004, Yang et al. 2006, Chan 
et al. 2008, Yeung et al. 2008, Vo et al. 2008, Furst et al. 2008, Garcfa-Lastra 
et al. 2010, Krasheninnikov et al. 2009). However, the difficulty is determining 
which specific molecules are present. In this study we show how changing the 
functionalization of the SWNT provides "another handle" for differentiating the 
SWNT's response to different gases/molecules. In this way, one may not only 
detect that a molecule is present, but also differentiate which molecule is present. 

Recent experimental advances now make the controlled doping of chirality 
selected SWNTs with metal atoms a possibility. Specifically, these include (1) 
photoluminesence, Raman and XAS techniques for measuring the fraction of 
various SWNT chiralities in an enriched sample (Kramberger et al. 2007, Ay- 
ala et al. 2009, De Blauwe et al. 2010); (2) the separation of SWNT samples 
by chirality using DNA wrapping (Zheng et al. 2003a, Zheng et al. 2003b, Tu 
et al. 2009, Li et al. 2007), chromatographic separation (Li et al. 2007, Fagan 
et al. 2007) and Density Gradient Ultracentrifugation (DGU) (Arnold et al. 2006); 
(3) SWNT resonators for measuring individual atoms of a metal vapor which ad- 
sorb on a SWNT (Lassagne et al. 2008, Jensen et al. 2008); and (4) aberration 
corrected low energy (<50 keV) transmission electron microscopy (TEM) (Jin 
et al. 2008). The latter provides control over the formation of defects in situ by 
adjusting the energy of the electron beam above and below the threshold energy 
for defect formation at a specific location on the SWNT. These methods provide 
such a high level of control that it is now possible for experimentalists to take a 
specific SWNT chirality and dope the structure with individual metal atoms at a 
specified location. 
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(6,6) SWNT Bandstructure 




Figure 3: Band energies s„^k relative to the Fermi level ep in eV versus momentum 
kinA~^ for a pristine (6,6) metallic armchair SWNT. 

At the same time, theorists are now able to embrace a "bottom up" approach 
to the design of nanosensors, harnessing the thermodynamics of self-assembly 
to find useful sensing devices in silico. With recent advances in both compu- 
tational power and methodologies, theorists can now efficiently and accurately 
screen hundreds of candidate sensor designs using a combination of DFT for en- 
ergetics of adsorption and stability, and NEGF methodologies for the electrical 
response (Garcfa-Lastra et al. 2010). 

In this chapter, we have used nine of the ten 3d TMs as candidates for substi- 
tutional doping of a (6,6) metallic armchair SWNT. The (6,6) SWNT was chosen 
as it is small enough to be efficiently calculated (24 atoms in the minimal unit 
ceh, 144 TM@MV or 143 atoms for the TM@DVI, TM@DVII), and at the same 
time experimentally realizable (Dresselhaus et al. 2001). We repeat the minimal 
unit cell 6 times to ensure convergence of the Hamiltonian to its bulk values at 
the borders of the cell. In other words, the unit cell is large enough that structural 
and electronic changes due to the TM in the vacancy are not felt at the unit cell 
boundary. The radius of the (6,6) SWNT is about 4. 1 A, while in experimental 
samples of chirality sorted tubes, the radius is typically around 7 A. 
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Carbon Nanotube Axis 
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Figure 4: Structural schematics and formation energies fiform in eV relative to a 
pristine (6,6) SWNT and a physisorbed transition metal atom for a 3d transition 
metal occupied monovacancy (TM@MV, black), divacancy I (TM@DVI, grey), 
and divacancy II (TM@DVII, white) in a (6,6) SWNT. Vacancy formation ener- 
gies for an empty monovacancy ( ), divacancy I ( ), and divacancy 

II ( ) are shown for comparison (Garcia-Lastra et al. 2010). 



However, in the -1 to 1 eV energy range of interest, for metallic tubes, the 
DOS is flat (2 eV~^), and has two open eigenchannels, so that T(e) = 2. This 
is clearly seen from the electronic band structure of the (6,6) SWNT, which is 
shown in Figure [3] In this way, the (6,6) SWNT results also provide a qualitative 
understanding of the response of larger metallic tubes. Further, for metallic arm- 
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chair (n, n) SWNTs both the transmission T(£) and DOS are completely flat near 
the Fermi level. This makes it much easier to differentiate the aff'ect of the TM 
dopant, and the adsorbed molecule on the transmission and DOS. 

The 3d TM atoms have been chosen as candidate dopants for several reasons. 
First, the 3d TM atoms are smaller for a given electronic configuration than the Ad 
and 5d TMs, and are expected to fit best in the structure of a SWNT, with the least 
strain. Even so, the TM atoms are typically pushed somewhat above the SWNT 
surface, due to both the strain on the graphene lattice, and the difference in atomic 
radii. These radii vary from the largest at 1.76 A for Ti, to the smallest at 1.42 
A for Ni (Ashcroft, and Mermin 1976). These are all between 2 and 3 times the 
atomic radius of C of 0.67 A. Also, the electronic properties of a TM dopant are 
mostly determined by the filling of the J-band, and should be quite similar for the 
Ad and 5d counterparts of the 3d TMs studied herein. 

The monovacancy, divacancy I and divacancy II were considered here since 
these are the three most stable vacancies, as depicted schematically in Figure |4j 
Comparing divacancy I and divacancy II results shows the effect of strain on the 
electronic properties of the active site, since for the divacancy II the strain is evenly 
distributed over the four bonds, while for the divacancy I the strain is mostly on 
two of the four bonds, as seen in Figure |4j For the monovacancy, the TM dopant 
is three fold coordinated, while for the divacancies, the TM dopant is four-fold 
coordinated. This explains why in general we find molecules are more strongly 
bonded to TMs in a monovacancy compared to either of the divacancies. 

For a strongly binding target molecule, such as CO, a divacancy provides the 
most suitable active site, as the site will not be completely saturated/covered, and 
a significant fraction of empty sites should always be present. On the other hand, 
for a weakly binding molecule, such as NH3, a monovacancy is more suitable, as 
there will then be a measurable change in coverage with molecule concentration. 



3 Density Functional Theory 

At its root, DFT is an exact reformulation of quantum mechanics in terms of the 
electron density (Kohn, and Sham 1965, Parr, and Yang 1989, Fiolhais et al. 2003, 
Engel, and Dreizler 2011). In this way it is a quantum mechanical method for 
calculating the ground state energy for many-body systems. From this theory, we 
may calculate formation energies Eads, adsorption energies Eads, and a system's 
Hamiltonian 'H and overlap S matrices, which form the first "puzzle piece" of 
our chemical nanosensor model, as shown in Figure |5j 
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^form[TM@VC] 

£„ds [X+TM@VC] 



5[Z+TM@VC] 
m^T] 5 [NT] 



DFT 





Figure 5: From density functional theory (DFT) we obtain the formation energies 
£^form of an active site (a 3d transition metal occupied vacancy (TM@VC)), the 
adsorption energies Eads for species X on the active site (X+TM@VC), and the 
Hamiltonian ^ and overlap matrices 5 for the scattering region with X adsorbed 
(X+TM@VC) and the leads (a pristine (6,6) SWNT (NT)). 

The reason why DFT is used is that it simplifies a problem in terms of the 
many electron wave function *P(ri, • • • ,r/v) into another problem, which is only 
in terms of the electron density p(ri). This is an important advantage because it 
reduces the coordinates from 3N to just 3 (x,y,z). 

¥(ri , . . . , Vn) 3N coordinates, 
eo[p] p(ri) 3 coordinates. 

This is based on the Hohenberg Kohn theorem which states that for a given exter- 
nal potential the ground state energy is a unique functional of the electron density. 




(3) 



where 



F[p] = T[p] + VM. 



(4) 
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Here T[p] is the kinetic energy, Vgelp] describes the electron-electron interaction 
and Vext is the external potential. However, an exact form for F[p] in terms of the 
density is unknown, so we approximate it by 

F[p] = T[p] + ^JJ jrJr'^^^ - J Jrp(r)v,,[p](r). (5) 

Here v^c is the exchange correlation potential that gives the correct electron-electron 
interaction Vee- v^c has corrections to the Coulomb potential in order to fulfill the 
Pauli exclusion principle, self interactions corrections, and electron-electron cor- 
relations. 

The most common implementation of DFT uses the Kohn-Sham method. The 
procedure for this method is to assume a form for Vxc, and then solve the non- 
interacting Schrodinger equation 

- — VV; + v(r)iA,- = Siif/i, (6) 
2m 

with the effective potential v given by 



v[p](r) = f dr'-f^ + vM + vM- (7) 
J l|r-r|l 

This potential depends on the density 

N/2 

p = J]ri^i, (8) 

/=0 

which in turn depends on the non-interacting Kohn-Sham wave functions i/r,, 
where is the number of electrons. We may thus obtain the wave functions and 
density by solving Eqns. (|6]), (|7]), and ([8]) in a self consistent manner. This Kohn- 
Sham self-consistency procedure is at the heart of DFT (Kohn, and Sham 1965). 

In this work, all total energy calculations and structural optimizations have 
been performed within the real- space density functional theory (DFT) code gpaw 
(Mortensen et al. 2005, Enkovaara et al. 2010) which is based on the projector 
augmented wave method (PAW) (Blochl 1994). This allows one to use a coarse 
grid, with "smoothed" pseudo-wave functions within the core, which describes 
the valence and conduction states quite well, to calculate the pseudo-density and 
pseudo-wave functions self-consistently. At the same time, it allows access to 
the full all-electron densities and wave functions, by projecting back onto the core 
wave functions on a finer grid. We find a grid spacing of 0.2 A provides a quite ac- 
curate representation of both the density and wave functions. For the purposes of 
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energetics of molecules and solids, the Perdew-Burke-Emzerhof implementation 
of the GGA for the exchange correlation functional (Perdew et al. 1996) works 
quite well. Spin polarization has been taken into account in all calculations. This 
is necessary, since 3d TM atoms often have a magnetic moment, as do molecules 
in the atmosphere such as O2. 

Pristine SWNTs are known to be chemically inert - a property closely related 
to their high stability. As a consequence, only radicals bind strong enough to 
the SWNT to notably affect its electrical properties (Charlier et al. 2007, Hierold 
2008, Valentini et al. 2004, ZanoUi, and Charlier 2009, Garcia-Lastra et al. 2010). 
To make SWNTs attractive for sensor applications thus requires some kind of 
functionalization, e.g. through doping or decoration of the SWNT sidewall (Pagan 
et al. 2003, Yagi et al. 2004, Yang et al. 2006, Ayala et al. 2007, Chan et al. 2008, 
Yeung et al. 2008, Vo et al. 2008, Fiirst et al. 2008, Garcia-Lastra et al. 2008, Ayala 
et al. 2008, Krasheninnikov et al. 2009, Ayala et al. 2010). Ideally, this type of 
functionalization could be used to control not only the reactivity of the SWNT but 
also the selectivity towards specific chemical species. 

Metallic doping of a (6,6) SWNT has been modeled in a supercell contain- 
ing six repeated minimal unit cells along the SWNT axis (dimensions: 15 Ax 15 
Ax 14.622 A). For this size of supercell a F-point sampling of the Brillouin zone 
was found to be sufficient to obtain reasonably converged adsorption and forma- 
tion energies. It should also be noted that for transport calculations within the 
NEGF method, it is implicitly assumed that the interaction between atoms near 
opposite ends of the unit cell are negligible, and any ^-point dependence along the 
transport direction is neglected. 

We define the formation energy for creating a vacancy (VC) occupied by a TM 
atom using the relation 

£form[TM@VC] = £[TM@VC] + nE[C] - E[TM + NT], (9) 

where £'[TM@VC] is the total energy of a TM atom occupying a vacancy in 
the nanotube, n is the number of carbon atoms removed to form the vacancy, 
^[C] is the energy per carbon atom in a pristine nanotube, and £'[TM-i-NT] is the 
total energy of the pristine nanotube with a physisorbed TM atom. The vacancies 
considered herein are the monovacancy and two divacancies shown in Figure |4j 
The energy required to form an empty vacancy is obtained from 

£form[C@VC] = E[YC] + nE[C] - E[NT], (10) 

where £'[VC] is the total energy of the nanotube with a vacancy of n atoms. 
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Table 1: Adsorption energies Ea^s in eV for N2, O2, H2O, CO, NH3, and H2S on a 
3d transition metal occupied mono vacancy (TM@MV). 





Ti 


V 


Cr 


Mn 


Fe 


Co 


Ni 


Cu 


Zn 


N2 


-0.61 


-0.76 


-0.57 


-0.70 


-0.77 


-0.60 


-0.49 


-0.34 


-0.51 


02 


-3.16 


-3.39 


-2.61 


-2.57 


-2.17 


-1.88 


-2.00 


-1.02 


-0.90 


H20 


-1.09 


-1.08 


-0.97 


-0.89 


-0.84 


-0.64 


-0.54 


-0.52 


-0.62 


CO 


-0.89 


-1.21 


-1.06 


-1.38 


-1.54 


-1.23 


-1.16 


-0.97 


-1.07 


NH3 


-1.39 


-1.46 


-1.35 


-1.32 


-1.31 


-1.03 


-0.94 


-0.89 


-1.12 


H2S 


-0.78 


-0.88 


-0.77 


-0.78 


-0.88 


-0.58 


-0.69 


-0.63 


-0.67 



Table 2: Adsorption energies E^^s in eV for N2, O2, H2O, CO, NH3, and H2S on a 
3d transition metal occupied di vacancy I (TM@DVI). 





Ti 


V 


Cr 


Mn 


Fe 


Co 


Ni 


Cu 


Zn 


N2 


-0.58 


-0.65 


-0.65 


-0.38 


-0.46 


-0.63 


-0.15 


-0.01 


-0.01 


02 


-2.72 


-3.41 


-2.74 


-1.92 


-1.67 


-1.42 


-0.61 


-0.06 


-0.12 


H20 


-1.21 


-1.19 


-0.79 


-0.49 


-0.51 


-0.33 


-0.19 


-0.05 


-0.25 


CO 


-0.80 


-1.11 


-1.14 


-1.06 


-1.62 


-1.90 


-1.10 


-0.12 


-0.04 


NH3 


-1.51 


-1.58 


-1.11 


-0.89 


-0.91 


-0.75 


-0.40 


-0.20 


-0.53 


H2S 


-0.87 


-0.78 


-2.13 


-0.50 


-0.64 


-0.59 


-0.31 


-0.04 


-0.17 



The calculated formation energies for the 3d TMs are shown in Figure|4j From 
the horizontal lines we see that both divacancies are more stable than the monova- 
cancy. This may be attributed to the presence of a two-fold coordinated C atom in 
the monovacancy, while all C atoms remain three-fold coordinated in the divacan- 
cies. When a TM atom occupies a vacancy, the strongest bonding to the C atoms 
is through its d orbitals (Griffith 1961). For this reason, Cu and Zn, which both 
have filled i^-bands, are rather unstable in the SWNT. For the remaining metals, 
adsorption in the monovacancies leads to quite stable structures. This is because 
the three-fold coordination of the C atoms and the SWNT's hexagonal structure 
are recovered when the metal atom is inserted. On the other hand, metal adsorp- 
tion in divacancies is slightly less stable because of the resulting pentagon defects, 
depicted in the upper panel in Figure |4} A similar behavior has been reported by 
Krasheninnikov et al. for TM atoms in graphene (Krasheninnikov et al. 2009). 

The adsorption energies for N2, O2, H2O, CO, NH3, and H2S on the metallic 
site of the doped (6,6) SWNTs are shown in Tables[T||2} and[3]for the TM occupied 
monovacancy, divacancy I, and divacancy II, respectively. The adsorption energy 
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Table 3: Adsorption energies Eads in eV for N2, O2, H2O, CO, NH3, and H2S on a 
3d transition metal occupied divacancy II (TM@DVII). 





Ti 


V 


Cr 


Mn 


Fe 


Co 


Ni 


Cu 


Zn 


N2 


-0.50 


-0.62 


-0.41 


-0.38 


-0.60 


-0.39 


-0.07 


-0.02 


-0.11 


02 


-2.27 


-2.64 


-2.35 


-2.17 


-1.48 


-1.31 


-0.53 


-0.09 


-0.18 


H20 


-1.15 


-1.20 


-0.74 


-0.64 


-0.49 


-0.55 


-0.16 


-0.05 


-0.17 


CO 


-0.78 


-1.10 


-1.05 


-0.96 


-1.49 


-1.30 


-0.82 


0.02 


0.01 


NH3 


-1.45 


-1.47 


-1.14 


-1.09 


-0.99 


-0.68 


-0.23 


-0.13 


-0.42 


H2S 


-0.80 


-0.83 


-0.62 


-0.62 


-0.56 


-0.50 


-0.10 


-0.01 


-0.10 



of a molecule X is defined by 

£ads[^+TM@VC] = E[X+TM@\C] - E[X] - £[TM@VC], (11) 

where £'[X@TM@ VC] is the total energy of molecule X on a TM atom occupying 
a vacancy, and E[X] is the gas phase energy of the molecule. 

From the adsorption energies shown in Tables [T| |2} and |3| we see that the 
earlier TMs tend to bind the adsorbates stronger than the late TMs. The latest 
metals in the series (Cu and Zn) bind adsorbates rather weakly in the divacancy 
structures. We also note that O2 binds significantly stronger than any of the three 
target molecules on Ti, V, Cr, and Mn (except for Cr in the divacancy I where H2S 
is found to spontaneously dissociate). Active sites containing these metals are 
therefore expected to be completely passivated if oxygen is present in the back- 
ground. Further, we find H2O is rather weakly bound to most of the active sites. 
This ensures that these types of sensors are robust against changes in humidity. 



4 Kinetic Modeling 

The next "piece of the puzzle" for our chemical nanosensor model is the descrip- 
tion of the coverage of our active sites by both target and background molecules 
using kinetic theory, as depicted in Figure [6} 

For a sensor containing hundreds, if not thousands of sites, the fractional cov- 
erage of sites is reasonably described by the fractional coverage in thermodynamic 
equilibrium at standard temperature and pressure in terms of the target molecule 
concentration. The adsorption reaction for a species X on an empty active site * 
is then 

X + *^X*. (12) 
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DFT 



Kinetics 



^form[TM@VC] 

£,ds [^+TM@VC] 

5[X+TM@VC] 
5 [NT] 




Figure 6: The energetics from density functional theory (DFT) are inputted to a 
kinetic model for the active site in thermodynamic equilibrium with concentra- 
tions Co, and gas phase entropies 5 gas, for each species X taken from experiment 
(Lide 2006-2007), to obtain the fractional coverage per active site (a 3d tran- 
sition metal occupied vacancy TM@VC) as a function of the target molecule's 
concentration C[X]. 



Under equihbrium conditions, the rate of adsorption and desorption of species X 
must be equal, so that 

radsm = kAX]@[*]C[X] - k.[X]@[X] = 0, (13) 

where rads is the rate of the adsorption reaction, k+ is the forward rate constant, 
k- is the backward rate constant, C is the atmospheric concentration, ©[*] is the 
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Table 4: Equilibrium atmospheric concentrations Co[X] at 4% humidity and gas 
phase entropies 5gas[^] in meV/K atT = 300 K where X is one of the background 
molecules N2, O2, and H2O or one of the target molecules CO, NH3, and H2S. 
Experimental gas phase entropies are taken from Ref. (Lide 2006-2007). 

X Com 5gasm X Com s,,,[x] 

N2 74.96% 1.988 me V/K CO 96.00 ppb 2.050 meV/K 
O2 20.11% 2.128 me V/K NH3 16.32 ppb 2.000 me V/K 
H2O 4.00% 1.959 me V/K H2S 0.96 ppb 2.136 meV/K 



fraction of active sites which are empty, and Q[X] is the fraction of active sites 
which are covered by species X. 



Using Eqn. ( 13 1 we may express the coverage of active sites by species X in 
terms of the ratio of forward and backward rate constants K = k+/k- and ©[*] as 



©m = i^m0[*]cm- 



(14) 



From kinetic theory, the ratio of rate constants may be obtained from the DFT 
adsorption energy Eads as 



K[X] = exp 



-Eads[X] - TS gas[X] 



(15) 



where ks is Boltzmann's constant, T is the temperature, and S gas[X] is the gas 
phase entropy of species X, as given in Table |4] (Lide 2006-2007). 

Since Eqn. ( [T3] ) is valid for all species in the set of background molecules S, 
so long as our description of the background is sufficient, if we sum up all the 
fractional coverages for each species, and the fraction of sites which are empty, 
we should obtain unity. More specifically. 



+ e[X] + 2]&[Y] = 1. 



(16) 



YeS 



By simple repeated substitution of Eqn. (14) into Eqn. (16 1, we may express 
the fractional coverage of active sites by species X in terms of the ratios of forward 
and backward rate constants, and concentrations in atmosphere as 

K[X]C[X] 



@[X] = 



(17) 



l+K[X]C[X] + ZYesK[Y]C[Y] 

To measure the "sensitivity" of an active site to the concentration of a target 
molecule, one should consider how the coverage of the active site by the target 
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Table 5: Logarithm of the change in fractional coverage with atmospheric concen- 
tration of the target molecule log [ffffj for CO, NH3, and H2S on a 3d transition 
metal occupied monovacancy (TM@MV). 





Ti 


V 


Cr 


Mn 


Fe 


Co 


Ni 


Cu 


Zn 


CO 


-86 


-82 


-58 


-44 


-22 


-23 


-30 


-7 


+9 


NH3 


-66 


-72 


-46 


-46 


-31 


-30 


-38 


-2 


+12 


H2S 


-91 


-96 


-70 


68 


-49 


-49 


-50 


-14 


-7 



Table 6: Logarithm of the change in fractional coverage with atmospheric concen- 
tration of the target molecule log [|§|fj], for CO, NH3, and H2S on a 3d transition 
metal occupied divacancy I (TM@DVI). 





Ti 


V 


Cr 


Mn 


Fe 


Co 


Ni 


Cu 


Zn 


CO 


-72 


-87 


-60 


-31 


+ 1 


+ 11 


+ 13 


-19 


-22 


NH3 


-44 


-68 


-60 


-37 


-27 


-28 


-10 


-15 


-2 


H2S 


-71 


-101 


-22 


-54 


-38 


-36 


-15 


-23 


-18 



Table 7: Logarithm of the change in fractional coverage with atmospheric concen- 
tration of the target molecule log [ff[fj], for CO, NH3, and H2S on a 3d transition 
metal occupied divacancy 11 (TM@DV11). 





Ti 


V 


Cr 


Mn 


Fe 


Co 


Ni 


Cu 


Zn 


CO 


-56 


-58 


-48 


-45 


+3 


+2 


+8 


-23 


-24 


NH3 


-29 


-43 


-44 


-39 


-16 


-21 


-14 


-18 


-7 


H2S 


-56 


-69 


-66 


-59 


-34 


-30 


-21 


-22 


-21 
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molecule changes with the concentration of the species in the atmosphere. It is 

only when there are changes in the active site's coverage, that through the result- 
ing changes in the sensing property, the nanosensor may be useful for a given 
species' detection. Specifically, we find the logarithm of the change in coverage 
with concentration of the target molecule, i.e. the derivative 



acts as a reasonable "descriptor" for an active site's sensitivity to species X. 
10° 
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Figure 7: Fractional coverage in thermal equilibrium of (a,b) Ni in a divacancy 
II (Ni@DVII) and (c,d) Cu in a monovacancy (Cu@MV) versus (a,c) CO concen- 
tration and (b,d) NH3 concentration in a background of air at room temperature 
and 1 bar of pressure (Mowbray et al. 2010). 
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As shown in Table [5| in the monovacancy, only the relatively inert Cu and 
Zn active sites are at all sensitive to the target molecules, specifically for CO and 
NH3. This is due to the binding energy of O2, which is > 2 eV for all other TM 
atoms occupying monovacancies, but only ~ 1 eV for both Cu and Zn, as shown 
in Table[T] This leaves "clean" Cu@MV and Zn@MV active sites, where CO or 
NH3 may adsorb. 

On the other hand, for both the divacancies shown in Tables |6] and |7] we find 
Fe, Co, and Ni are sensitive to CO. We attribute this to the similarity in binding 
energies for O2 and CO on these sites, which are between 0.5 and 1 .6 eV, as shown 
in Tables |2] and |3j We again find Ti, V, Cr, and Mn are completely oxidized, as 
was the case in the monovacancy, making these sites inactive, while for Cu and 
Zn the binding energies of the target molecules are too weak for them to adsorb 
on these sites. 

In Figure |7] we see the coverages of both Ni and Cu active sites as a function of 
CO and NH3 concentration. Clearly, both Cu and Ni active sites will be sensitive 
to the adsorption of CO. On the other hand, the coverage of the Cu active sites is 
sensitive to NH3 concentration, while Ni active sites are not. This suggests that by 
combining the response of Cu and Ni active sites, we may obtain a multifunctional 
sensor for both CO and NH3. 

5 Non-Equilibrium Green's Function Methodology 

The third "puzzle piece" needed to model our chemical nanosensor is the descrip- 
tion of the sensing property, specifically the change in resistance of the sensor 
using the NEGF method, as depicted in Figure [8| Although molecules may ad- 
sorb on the active sites, the sensing property of the SWNT resistance must change 
significantly for the sensor to be useful. In the following sections we show both 
the density of states (DOS) and transmission probability T(e) for an electron of a 
given energy e to travel past an active site as calculated within the NEGF method. 

Since the energy levels of a given molecule may be thought of as its elec- 
tronic "finger print", so long as there is sufficient binding between the molecule 
and active site, this "finger print" is evidenced in the DOS of the system. Further- 
more, the peaks of the DOS due to the molecule tend to form Fano anti-resonances 
(Fiirst et al. 2008) in the transmission, as electrons are scattered off these states. 
In this way, any adsorbed molecule tends to leave its "finger print" in the trans- 
mission probability through an active site. For this reason, conductance/resistance 
is generally an effective sensing property. 
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DFT 



^form[TM@VC] 

£,ds [^+TM@VC] 

5[X+TM@VC] 
5 [NT] 



NEGF 



T[Z+TM@VC] (8) 
DOS[Z+TM@VC] (e) 



Kinetics 





J 



Figure 8: The Hamiltonian and overlap 5 matrices obtained from density func- 
tional theory (DFT) are used within the non-equilibrium Green's function (NEGF) 
methodology to obtain the density of states (DOS) and transmission probabihty T 
for each species X on an active site (X-i-TM@VC), as a function of energy s. 



The NEGF methodology (Datta 1997) is based on the Green's function, which 
is dependent on two space-time coordinates. The NEGF method can be applied to 
both extended and infinite systems and it can handle strong external fields without 
being perturbed. The approximations made in the NEGF method can be cho- 
sen in order to satisfy macroscopic conservation laws. An interesting feature of 
the NEGF method is that memory effects in transport that occur due to electron- 
electron interactions can be analyzed. 

As depicted schematically in Figure [8} the NEGF method requires the Hamil- 
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Ionian and overlap S matrices for both the clean and doped SWNT struc- 
tures. The NEGF calculation is then performed using as "leads" the clean SWNT 
Hamiltonian, with a "scattering region" defined by the doped SWNT Hamiltonian. 
Within this representation, the states which form the Hamiltonian must be in some 
way localized, so that the Hamiltonian of the scattering region is "converged" to 
the Hamiltonian of the leads near the unit cell boundaries. This may be accom- 
plished through the use of Wannier functions to localize periodic wave functions, 
or a separate calculation with a locally centered atomic orbital (LCAO) basis set 
(Strange et al. 2008). 

In this work, transport calculations for the optimized structures have been 
performed with an electronic Hamiltonian obtained from the siesta code (Soler 
et al. 2002) in a double zeta polarized (DZP) LCAO basis set. Previous studies 
have shown that a DZP basis set is sufficient to converge the Hamiltonian, yield- 
ing similar results to those obtained from a plane- wave calculation using Wannier 
functions (Strange et al. 2008). 

A large supercell was thus necessary for the Hamiltonian, so that each of the 
four SWNT layers adjacent to the boundaries were within 0.1 eV of the 

Hamiltonian for the respective leads :?4, i-e. max|i^"" - < 0.1 eV. In this 
way the electronic structure at the edges of the central region was ensured to be 
converged to that in the leads. 

The Landauer-Biitticker conductance for a multi-terminal system can be calcu- 
lated from the Green's function of the central region, Gc, according to the formula 
(Meir, and Wingreen 1992, Thygesen 2006, Strange et al. 2008) 



where the trace runs over all localized basis functions in the central region, and 
Go = 2e^/h is the quantum of conductance. To describe the conductance at small 
bias for semiconducting systems, the Fermi energy sp should be taken as the en- 
ergy of the valence-band maximum svb or conduction-band minimum scb for p- 
type and n-type semiconductors, respectively. The central region Green's function 
is calculated from 



where z = s + iO'^, Sc and ^Hc are the overlap matrix and Kohn-Sham Hamiltonian 
matrix of the central region in the localized basis, Ha is the self-energy of lead a. 



G = GoTr[G, 



C^inGf-'Tout. 



(19) 




(20) 



(21) 
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and the coupling elements between the central region and lead a for the overlap 
and Kohn-Sham Hamiltonian are Sca and ^Hca respectively. 

The coupling strengths of the input and output leads are then given by Ym/out = 

In the following sections, we will analyze the transmission data for the clean 
TM doped SWNTs. To aid our understanding we will also consider the DOS of 
the system, which drives the behavior we see in the transmission. We will first 
analyze all the different metals in the divancancy II, then the same will be done 
with the divacancy I and the monovacancy. Since we will concentrate on the DOS, 
it is important to understand the way the DOS and transmission T(e) are related. 
The peaks found in the DOS match with "dips" in the transmission, which are 
Fano anti -resonances (Fiirst et al. 2008). While narrow peaks in the DOS denote 
localized states, broad peaks denote strong bonding or hybridized states. 



5.1 Divacancy II 

For a TM occupying a divacancy II (DVII) in a (6,6) SWNT (TM@DVII), shown 
schematically in Figure |4| there is an octahedral bonding of the TM atom, similar 
to the bonding in rutile metal oxides. Here, the TM bonding with the SWNT is 
through the two hybridized 3d orbitals of the TM. From Figure |4j we see 

that all the TM-C bonds are equally strained by the curvature of the nanotube. 
This means that all TM-C bonds are equivalent. To put it in other words, there is 
one type of four-fold degenerate TM-C bonding state. 

The DOS and transmission T(e) through a pure TM occupied divacancy (TM@ 
DVII) are shown in Figures [9] and [10} respectively. As mentioned, if we compare 
the DOS and T(e) shown in Figures[9]and[10[ we find peaks and dips, respectively, 
relative to the flat DOS and T(e) of a pristine (6,6) SWNT. 

For Ti@DVII there are four occupied C n — Ti "^''^'" bonds. These bonds are 
more than 1 eV below the Fermi level, so that they are not relevant for this study. 
There are also present completely unoccupied (i^.2_j,2 orbitals and anti-bonding Ti 
-^^^ orbitals. However, only the d orbitals of the TMs are important in the -1 
to 1 eV range. These d orbitals are what determines the binding and transmission 
for the TM@DVII system. 

In the case of V@DVII, there is one extra electron in the d^2_y2 orbital, which 
is a non-bonding orbital. This extra electron yields a spin splitting in the DOS, 
as seen in Figure |9| The d^2_y2 orbital is, by symmetry, non-bonding and it is 
localized on the V atom. 
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(c) Cr@DVII 




I I I I I I I I 



(f) Co@ DVII 



(i) Zn@DVII 



J I I I I I I I I 

1 1 



8 - £p (eV) 



Figure 9: Total density of states (DOS) in eV ' versus energy e in eV relative to 
the Fermi level £/r for a 3J transition metal (TM) occupied divacancy II (DVII) in 

a (6,6) metallic armchair SWNT (TM@DVII). Results for spin T ( ), spin i (- 

— ), and spin averaged ^ ( ) are shown. The constant DOS at 2 eV"' for a 

pristine (6,6) SWNT is provided for reference. 



For Cr@DVII there are two extra electrons in the d^i.f. orbital, instead of only 
one. The states are shifted down in energy and there is no state around the Fermi 
level, as seen in Figure |9} 

With Mn@DVII, electrons begin occupying the anti-bonding orbitals. There 
are then localized spin states above the Fermi level. These spin unpaired d states, 
which are localized on the TM atom, move down in energy as they are filled with 
more J-electrons. This is what we see in Figure|9]for Fe@DVII and Co@DVII. 

For Ni@DVII, these localized states become spin-paired, with one of them 
being completely occupied, and the other one completely empty. When we add 
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Figure 10: Transmission versus energy s in eV relative to the Fermi level sp 
through a 3d transition metal (TM) occupied divacancy II (DVII) in a (6,6) metal- 
lic armchair SWNT (TM@DVII). Results for spin T (• • ■ ■), spin i ( ), and 

spin averaged ^ ( ) are shown. The constant transmission of 2 for a pristine 

(6,6) SWNT is provided for reference. 



another electron to the d band, moving to Cu@DVII, both states become fully 
occupied, localized on the Cu atom, and spin-paired. 

Finally, for Zn@DVII, the J-orbitals are now completely filled. The DOS 
closely resembles that of Ti@DVII, and the states are spin-paired, as shown in 
Figure |9j 



26 



10 



> 

CD 

o 

Q 



> 

CD 

O 
Q 



> 
CO 

o 



(a) Ti@DVI 





_] I I I I I I I L. 





£-£p (eV) 



(b) 




\/@DVI 




J I I I I I I I L_ 



(h) Cu@DVI 



_] I I I I I I I L. 



1 -1 





£ - £p (eV) 





_l I I I I I I I L 



(i) Zn@DVI 




_] I I I I I I I L. 



£ - £p (eV) 



Figure 11: Total density of states (DOS) in eV ^ versus energy s in eV relative to 
the Fermi level Sf for a 3d transition metal (TM) occupied divacancy I (DVI) in 

a (6,6) metallic armchair SWNT (TM@DVI). Results for spin T ( ), spin i (- 

— ), and spin averaged ^^ ( ) are shown. The constant DOS at 2 eV~' for a 

pristine (6,6) SWNT is provided for reference. 



5.2 Divacancy I 

For a TM occupying a divacancy I (DVI) in a (6,6) SWNT (TM@DVI), shown 
schematically in Figure]?} the states are no longer as localized on the metal atom 
as for the TM@DVII, and are more strongly hybridized with the C n states. Due 
to the geometric configuration of TM@DVI, the angle of two of the bonds with 
the axis of the nanotube is much smaller than the angle of the other two bonds. 
We will call these two different bonds near-axis and off"-axis, respectively. This 
difference in the orientation produces a difference in the strain the bonds will have 
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Figure 12: Transmission versus energy s in eV relative to the Fermi level sp 
through a 3d transition metal (TM) occupied divacancy I (DVI) in a (6,6) metallic 

armchair SWNT (TM@DVI). Results for spin T (■ • • 0, spin i ( ), and spin 

averaged ^ ( ) are shown. The constant transmission of 2 for a pristine (6,6) 

SWNT is provided for reference. 



to withstand. 

The off-axis bonds will undergo more stress due the SWNT's curvature, com- 
pared to the near-axis bonds. This means that there are two types of degenerate 
TM-C bonding states. This degeneration results in a larger number of broader 
peaks and dips, compared to those obtained for a TM@DVII, as shown in Figures 



1 1 and [121 respectively 



Even with these differences it is still possible to see in Figure 11 the filling 
of states, and their movement to lower energy, exactly like in the case of the 
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Figure 13: Total density of states (DOS) in eV ^ versus energy s in eV relative to 
the Fermi level Sf for a 3d transition metal (TM) occupied monovacancy (MV) in 

a (6,6) metallic armchair SWNT (TM@MV). Results for spin T ( ), spin i (- 

— ), and spin averaged ^^ ( ) are shown. The constant DOS at 2 eV"' for a 

pristine (6,6) SWNT is provided for reference. 



TM@DVII. It is very likely that the differences in the strain contribute to the 
splitting of the spin degeneracy for Ti@DVI and Ni@DVI. 

The Ti@DVI d orbitals are half-filled through bonding with the C n orbitals, 
as is the 4^ level. The Cu@DVI d orbitals are filled by taking an electron from the 
45 level. So, in the case of Cu@DVI, bonding occurs through the half-occupied 
4s orbital. All of the Zn d orbitals are completely filled, as is the 4^' level. For 
these three systems, the DOS shown in Figure [TT] is very similar to that shown 
in Figure |9] for the TM@DVII case. However, for the Ti@DVI DOS there is a 
noticeable difference, most probably due to the differences in strain for the two 
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Figure 14: Transmission versus energy e in eV relative to the Fermi level sp 
through a 3d transition metal (TM) occupied monovacancy (MV) in a (6,6) metal- 
lic armchair SWNT (TM@MV). Results for spin T ( ), spin [{ ), and spin 

averaged ^ ( ) are shown. The constant transmission of 2 for a pristine (6,6) 

SWNT is provided for reference. 



types of C-TM bonds (near- axis and ofF-axis). 



5.3 Monovacancy 

For a TM occupying a monovacancy (MV) in a (6,6) SWNT (TM@MV), shown 
schematically in Figure |4} the TM atom has three bonds with three C atoms, two 
of the bonds being shorter, stronger, and near-axis, and one longer, weaker, and 
perpendicular to the nanotube axis. Not only is this bond weaker, but it is also 
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Figure 15: Total density of states (DOS) in eV ' versus energy s in eV relative to 
the Fermi level sp for a 3d transition metal (TM) occupied monovacancy (MV) 
in a (6,6) metallic armchair SWNT (TM@MV). Results for the clean TM@MV 

( ) and with adsorbed N2 ( ), O2 ( ), H2O ( ), CO ( 

-), and NH3 (- • — •) are shown. The constant DOS at 2 eV"' for a pristine (6,6) 
SWNT is provided for reference. 



the one that withstands the largest strain. From Figure 13 we see that for Ti@MV 
and Ni@MV we have spin splitting, which may be attributable to the strain in the 
weaker C-TM bond. Figure [13] also shows that for the monovacancy the DOS 
shifts down in energy as the d orbitals on the TM atom are filled. The differences 
between the bond lengths and strengths of the three C-TM bonds yields several 
broad peaks in the DOS, as seen in Figure [131 and likewise broad dips in the 



transmission, as seen in Figure 14 



As shown in Figure n3l for Ti@MV, V@MV, and Cr@MV, we see a broad 
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Figure 16: Transmission versus energy s in eV relative to the Fermi level sp 
through a 3d transition metal (TM) occupied monovacancy (MV) in a (6,6) metal- 
lic armchair SWNT (TM@MV). Results for the clean TM@MV ( ) and 

with adsorbed N2 ( ), O2 ( ), H2O ( ), CO ( ), and NH3 

(-•■-■■) are shown. The constant transmission of 2 for a pristine (6,6) SWNT is 
provided for reference. 



spin-unpaired peak centered near the Fermi level. This peak is shifted down in 
energy for Mn@MV, which has five d electrons. At this point a narrower peak, 
for a J-state more localized on the TM atom enters the energy window and it is 
positioned at round 0.2 eV. This state broadens for Co@MV where it is partly oc- 
cupied, and nearly spin-paired. For Ni@MV this state goes further down in energy 
relative to the Fermi level and becomes spin paired for Cu@MV and Zn@MV. 
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Figure 17: Total density of states (DOS) in eV ^ versus energy s in eV relative 
to the Fermi level sp for a 3d transition metal (TM) occupied divacancy I (DVI) 
in a (6,6) metallic armchair SWNT (TM@DVI). Results for the clean TM@MV 

( ) and with adsorbed N2 ( ), O2 ( ), H2O ( ), CO ( 

-), and NH3 (- • — •) are shown. The constant DOS at 2 eV"' for a pristine (6,6) 
SWNT is provided for reference. 



5.4 Target and Background Molecules 

To estimate the effect of adsorbates on the electrical conductance of doped SWNTs, 
we first consider the change in conductance when a single molecule is adsorbed 
on a metal site of an otherwise pristine SWNT. In Figs. T5| 17} and 19 we show the 
calculated DOS for a TM occupied monovacancy, divacancy I, and divacancy II, 
respectively, with and without an adsorbed molecule X. These may be compared 
with Figs. 16} 18, and 20 where we show the calculated transmission probability 
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Figure 18: Transmission versus energy s in eV relative to the Fermi level sp 
through a 3d transition metal (TM) occupied divacancy I (DVI) in a (6,6) metallic 

armchair SWNT (TM@DVI). Results for the clean TM@DVI ( ) and with 

adsorbed N2 ( ), O2 ( ), H2O ( ), CO ( ), and NH3 (- • 

• - • •) are shown. The constant transmission of 2 for a pristine (6,6) SWNT is 
provided for reference. 



T for a TM occupied monovacancy, divacancy I, and divacancy II, respectively, 
with and without an adsorbed molecule X. 

In contrast to the adsorption energies, there are no clear trends in the conduc- 
tances. The sensitivity of the conductance is perhaps most clearly demonstrated 
by the absence of correlation between different types of vacancies. 

Close to the Fermi level, the conductance of a perfect armchair SWNT equals 
2Go. The presence of the metal dopant leads to several dips in the transmission 
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Figure 19: Total density of states (DOS) in eV ^ versus energy e in eV relative 
to the Fermi level sp for a 3J transition metal (TM) occupied divacancy II (DVII) 
in a (6,6) metallic armchair SWNT (TM@DVII). Resuks for the clean TM@MV 

( ) and with adsorbed N2 ( ), O2 ( ), H2O ( ), CO ( 

-), and NH3 (- • — •) are shown. The constant DOS at 2 eV"' for a pristine (6,6) 
SWNT is provided for reference. 



function known as Fano anti-resonances (Fiirst et al. 2008). The position and 
shape of these dips depends on the J-levels of the TM atom, the character of its 
bonding to the SWNT, and is further affected by the presence of the adsorbate 
molecule. The coupling of all these factors is very complex and makes it difficult 
to estimate or rationalize the value of the conductance. For the spin-polarized 
cases, we use the spin-averaged conductances, i.e. G = (G^ -v G\)I2. 

Cu@DVI binds all molecules rather weakly, so any changes in the DOS from 
the clean Cu@DVI structure between -1 and 1 eV are attributable to states com- 
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Figure 20: Transmission versus energy s in eV relative to the Fermi level sp 
through a 3d transition metal (TM) occupied divacancy II (DVII) in a (6,6) metal- 
lic armchair SWNT (TM@DVII). Results for the clean TM@DVII ( ) and 

with adsorbed N2 ( ), O2 ( ), H2O ( ), CO ( ), and NH3 

(-•■-■■) are shown. The constant transmission of 2 for a pristine (6,6) SWNT is 
provided for reference. 



pletely localized on the molecule. Here we only see two states localized on the 
O2 molecule between 0.6 and 0.8 eV above the Fermi level. For CO-i-Cu@DVI 
we see a shifting of the DOS, and hybridization with the Cu@DVI states. In par- 
ticular, the broad peak below Sf in the Cu@DVI DOS, which is binding with the 
C n states, is shifted down in energy by about 0. 1 eV when CO is adsorbed on the 
Cu@DVI site. 

For O2, its ground state has spin 2, with one electron occupying a — anti- 
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bonding n* orbital, and the other electron occupying a Py — Py anti-bonding n* 
orbital (Mowbray et al. 2009). For O2 on Ti@DVI, V@DVI, and Cr@DVI, the 
Py — p, anti-bonding n* orbital hybridizes with the transition metal atom's dy, state, 
leaving the spin down component of the other py — Py anti-bonding n* state on the 
molecule unoccupied. This state is visible as a sharp peak above the Fermi level. 
By symmetry the Py — Py anti -bonding O2 n* state cannot bond to the TM@DVI 
system. 

It should be noted that O2 is the only molecule with molecular states in the 
-1 to -1-1 eV range relative to the (6,6) SWNT Fermi level. The O2 Py — Py anti- 
bonding n* states by symmetry are forbidden from bonding with the TM@DVI 
and TM@DVII systems. On the other hand, the O2 Pz — Pz anti-bonding n* state 
overlaps with the TM dyz state. For Cu and Zn, as the d orbitals are filled, we find 
the O2 anti-bonding n* states are both spin unpaired, and exhibit little hybridiza- 
tion with the TM@VC site. 

Overall, the TM d states are shifted up and down in energy when a molecule 
is adsorbed, based on both charge transfer, and the symmetry of the molecule's 
orbitals. This leads to a rich response in the DOS of the TM@VC system upon 
adsorption of a molecule. The main conclusion to be made is that the DOS and 
thus the transmission through the TM@VC systems are quite sensitive to the ad- 
sorption of a molecule, and give a strong response. This makes the resistance 
through such active sites an excellent sensing property. 



6 Sensing Property 

The final "piece of the puzzle" for our model of a chemical nanosensor is the 
description of the sensing property, i.e. the change in resistance upon adsorption 
of a molecule under atmospheric conditions, via basic circuit theory, as depicted 



in Figure 21 At this point we bring together the transmission T(e) calculated from 
the NEGF method and the coverage of active sites &[X] from the kinetic model 
to describe the change in resistance AR through the doped SWNT as a function of 
the target molecule's concentration C[X]. 

We now estimate the resistance of a SWNT containing several impurities (a 
specific metal dopant with different molecular adsorbates). Under the assumption 
that the electron phase-coherence length, -d^, is smaller than the average distance 
between the dopants, d, we may neglect quantum interference and obtain the total 
resistance by adding the scattering resistances due to each impurity separately. 
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Table 8: Conductance G[X] in Go = 2e^/h for a 3d transition metal occupied 
monovacancy (TM@MV) with N2, O2, H2O, CO, NH3, and H2S adsorbed. The 
conductance of the clean TM@MV is provided for comparison. 





Ti 


V 


Cr 


Mn 


Fe 


Co 


Ni 


Cu 


Zn 


N2 


1.44 


1.63 


1.44 


1.72 


1.88 


1.19 


1.40 


0.97 


1.16 


02 


1.44 


1.78 


1.72 


1.66 


1.56 


1.83 


1.33 


1.13 


1.69 


H20 


1.42 


1.31 


1.64 


1.39 


1.50 


1.47 


1.18 


0.97 


0.99 


CO 


1.31 


1.61 


1.61 


1.59 


1.33 


1.09 


1.58 


0.97 


1.79 


NH3 


1.57 


1.37 


1.68 


1.47 


1.51 


1.50 


1.39 


1.05 


1.01 


H2S 


1.40 


1.51 


1.64 


1.51 


1.50 


1.43 


1.32 


0.98 


1.13 


Clean 


1.40 


1.41 


1.38 


1.77 


1.87 


1.08 


1.06 


1.31 


1.02 



Table 9: Conductance G[X] in Go = 2e^/h for a 3d transition metal occupied 
divacancy I (TM@DVI) with N2, O2, H2O, CO, NH3, and H2S adsorbed. The 
conductance of the clean TM@DVI is provided for comparison. 





Ti 


V 


Cr 


Mn 


Fe 


Co 


Ni 


Cu 


Zn 


N2 


1.52 


1.69 


1.60 


1.38 


1.20 


1.16 


1.06 


1.28 


1.42 


02 


1.13 


1.24 


1.54 


1.51 


1.48 


1.09 


1.33 


1.28 


1.30 


H20 


1.50 


1.52 


1.69 


1.27 


1.49 


1.07 


1.83 


1.25 


1.44 


CO 


1.58 


1.29 


1.28 


1.50 


1.74 


1.13 


1.56 


1.50 


1.48 


NH3 


1.50 


1.25 


1.68 


1.27 


1.42 


1.14 


1.05 


1.35 


1.51 


H2S 


1.51 


1.59 


1.69 


1.31 


1.43 


1.11 


1.03 


1.29 


1.47 


Clean 


1.61 


1.65 


1.73 


1.29 


1.57 


1.02 


1.57 


1.25 


1.34 



Table 10: Conductance G[X] in Go = 2e^/h for a 3d transition metal occupied 
divacancy II (TM@DVII) with N2, O2, H2O, CO, NH3, and H2S adsorbed. The 
conductance of the clean TM@DVII is provided for comparison. 





Ti 


V 


Cr 


Mn 


Fe 


Co 


Ni 


Cu 


Zn 


N2 


1.65 


1.77 


1.65 


1.27 


1.12 


1.33 


0.51 


1.06 


1.49 


02 


1.65 


1.29 


1.47 


1.54 


1.15 


1.49 


1.00 


1.04 


1.34 


H20 


1.57 


1.67 


1.69 


1.20 


1.21 


1.45 


1.29 


1.02 


1.56 


CO 


1.68 


1.79 


1.52 


1.29 


1.07 


0.61 


0.69 


1.17 


1.56 


NH3 


1.61 


1.69 


1.67 


1.19 


1.18 


1.50 


0.86 


1.01 


1.60 


H2S 


1.62 


1.71 


1.65 


1.23 


1.15 


1.43 


0.79 


1.06 


1.58 


Clean 


1.60 


1.78 


1.71 


1.22 


1.06 


1.41 


1.71 


1.04 


1.49 



38 



DFT 



Kinetics 



^form[TM@VC] 

£,ds [^+TM@VC] 

5[X+TM@VC] 
5 [NT] 




T[X+TM@VC] (£)N^ 
DOS[Z+TM@VC] (e) 




Figure 21: The coverage and transmission probabilities T(e) are used within 
simple circuit theory to estimate the scattering resistance Rs for species X on an 
active site (X+TM@VC), and hence derive the overall change in resistance AR 
for a particular active site (TM@VC) as a function of the target molecule's con- 
centration C[X]. 



The scattering resistance due to a single impurity is then given by 

R^[X] = — —. (22) 

G[X] 2Go 

Here G[X] = r[X-l-TM@VC](eF) is the Landauer conductance of the pristine 
SWNT with a single metal dopant occupied by molecule X, i.e. the transmission 
at the Fermi level sp, l/(2Go) is the contact resistance of a (6,6) SWNT, and 
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Figure 22: Change in resistance AR in (a,b) Q and (c,d) per active site versus 
(a,c) CO and (b,d) NH3 concentration in ppm for Ni in a divacancy II (Ni@DVII,- 

) and Cu in a monovacancy (Cu@MV, ) of a (6,6) carbon nanotube. 

The reference concentrations of CO and NH3 are 0.1 and 0.01 ppm, respectively, 
in a background of air at room temperature and 1 bar of pressure. 



Go = 2e^/h is the quantum of conductance. The calculated conductances are 
provided in Tables [8| [9} and [10] for a molecule X adsorbed on a 3d TM occupied 



monovacancy, divacancy I, and divacancy II, respectively. 

We may now obtain the average change in resistance AR of an active site as 
a function of target molecule concentration. As discussed in Ref. (Garcfa-Lastra 
et al. 2010), this change in resistance is reasonably well described by 

AR^Yj ^.vm(0K C] - 0[X, Co]), (23) 

X 
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where C is the concentration, and Q is the concentration at standard temperature 
and pressure, as given in Table |4} 

In Figure 22 we show AR for a single active site (Ni or Cu) as a function of 
target molecule concentration (CO or NH3). Keeping in mind that a change in 
concentration from 1 ppm to 50 ppm amounts to a change from allowable to toxic 
concentrations for both CO and NH3, we see that this sensor design should work 
effectively for both target molecules. 



7 Conclusions 

In summary, we have presented a general model for the design of chemical nanosen- 
sors which takes the adsorption energies of the relevant chemical species and their 
individual scattering resistances as the only input. On the basis of this model we 
have performed a computational screening of TM doped SWNTs, and demon- 
strated using ab initio calculations how a combined Cu and Ni doped metallic 
SWNT device may work effectively as a multifunctional sensor for both CO and 
NH3. Further, we find that by varying the metal dopant, we obtain "another han- 
dle" for tuning the sensitivity of our sensor. The methodology employed may also 
be applied to other nanosensor designs and different environments, and demon- 
strates the potential of computational screening studies for the design of chemical 
nanosensors in silico. 
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